Effects of trap anisotropy on impurity scattering regime in a Fermi gas 
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We evaluate the low-lying oscillation modes and the ballistic expansion properties of a harmon- 
ically trapped gas of fermionic 40 K atoms containing thermal 87 Rb impurities as functions of the 
anisotropy of the trap. Numerical results are obtained by solving the Vlasov-Landau equations 
for the one-body phase-space distribution functions and are used to test simple scaling Ansatzes. 
Starting from the gas in a weak impurity-scattering regime inside a spherical trap, the time scales 
associated to motions in the axial and azimuthal directions enter into competition as the trap is 
deformed to an elongated cigar-like shape. This competition gives rise to coexistence of collisionless 
and hydrodynamic behaviors in the low-lying surface modes of the gas as well as to a dependence 
of the aspect ratio of the expanding cloud on the collision time. 
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I. INTRODUCTION 

Experiments on ballistic expansion and collective- 
mode excitations have provided important diagnostic 
tools in the study of quantum-degenerate atomic and 
molecular gases (for a recent review see Q). A mile- 
stone example is the interpretation of the anisotropic ex- 
pansion of a cloud of ultracold bosons as the smoking 
gun of Bose-Einstein condensation in a trapped dilute 
gas 0JS 0] • The later achievement of degenerate Fermi 
gases U and boson- fermion mixtures [flUEl nas opened 
up the possibility to experimentally observe rich quan- 
tum phase diagrams and to study the interplay between 
external confinement and interspecies collisions. More 
recently, the expansion behavior of strongly interacting 
fermion clouds has been the object of experiments aimed 
at probing novel superfluid states on the crossover from 
a BCS superfluid to a condensate of molecular dimers 

II US El. 

Collective modes provide a direct measure of the col- 
lisionality of a dilute quantum gas. Monopolar and 
quadrupolar modes have been proposed as markers of the 
approach of quantum phase transitions in binary mix- 
tures of Bose-Einstein condensed gases (BEC's) |l2|, in 
boson- fermion mixtures |T3.ll4j, and in fermion mixtures 
across the BCS-BEC crossover [HE [HI. The dynami- 
cal transition from collisionless to hydrodynamic behav- 
ior in fermion mixtures has been followed experimentally 
|18| and numerically pjjj in two-component fermion mix- 
tures by studying their dipolar oscillation modes. The 
role played in this context by mobile impurities inside a 
fermion gas under spherical confinement has also been 
studied by numerical means po| . 

In the present work we examine how the anisotropy 
of the trap affects the low-lying oscillation modes and 
the ballistic expansion of a gas of fermionic 40 K atoms 
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containing a small concentration of thermal Rb atoms. 
We numerically solve the Vlasov-Landau equations for 
the evolution of the phase-space distribution functions 
within a particle- in-cell approach and compare the results 
with simple scaling Ansatzes. Our results demonstrate 
that in a cigar-shaped harmonic confinement, where two 
different confinement frequencies determine two different 
time scales, collisionless and hydrodynamic behaviors can 
coexist in the low-lying collective excitations of the gas 
and that the aspect ratio of the expanding cloud shows a 
non-monotonic dependence on the anisotropy of the trap. 

The paper is organized as follows. In Sect. IJwe in- 
troduce the mixture under study and the basic equations 
that will be used to describe its dynamics. Section IIIII 
analyzes the monopolar, dipolar, and quadrupolar oscil- 
lations of the gas as functions of the anisotropy of the 
confining potential, while Sect. II VI gives a discussion of 
the dynamics of the free expansion of the fermion cloud. 
Finally, Sect. presents a summary and the main con- 
clusions of our work. 



II. THE MIXTURE 

We consider a gas of fermionic atoms of mass m p con- 
fined inside an axially symmetric harmonic trap of the 
form 

V F (r) = ^m F LU F ^(x 2 + y 2 + \ 2 z 2 ), (1) 

where u>p± is the angular trap frequency for motions 
along the x and y directions and A = LUp^ z /ujp± is the 
anisotropy parameter. Potentials with A < 1 generate 
cigar-shaped density profiles. The fermions are wholly 
spin-polarized and hence at very low temperature the 
Pauli principle quenches collisions among them. A sec- 
ond component must then be added in order to thermal- 
ize the gas and to drive its collisionality. 

We dope the fermion gas with a small number of bosons 
having larger atomic mass m b ■ We simply think of them 
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as impurities confined by the harmonic potential 

V B (r)^^m B uj 2 B ^(x 2 + y 2 + X 2 z 2 ), (2) 

where u>b,± is the angular trap frequency in the az- 
imuthal plane and we have assumed that bosons and 
fcrmions share the same trap anisotropy A. The dif- 
ference in atomic masses and particle numbers causes 
a much lower quantum degeneracy temperature for the 
bosons and we can therefore consider them as an uncon- 
densed cloud, even in the presence of highly degenerate 
fermions. 

The dynamics of the mixture is described through 
the one-body distribution functions f^ F ' B '(r,p,t) in the 
Boltzmann approximation. Their evolution is governed 
by the Vlasov-Landau kinetic equations (VLE), 

d t f U) + — •V r /«-V r C/«-Vp/« =C[f {F \f {B ^ (3) 

where the Hartree-Fock effective potential is U^> (r, t) — 
Vj(r) + gn^\r,t) with j denoting the species different 
from j. Here we have set g = 2irh a/m r with a being the 
s-wave scattering length of a fermion-boson pair and m r 
its reduced mass, and rft' (r, t) is the spatial density given 
by integration of f^'(r,p,t) over the momentum de- 
grees of freedom. Since we deal with low concentrations 
of impurities we have neglected impurity-impurity inter- 
actions. In addition, collisions between spin-polarized 
fermions are negligible at low temperature and thus the 
collision integral C in Eq. J3J involves only collisions be- 
tween fermions and impurities. This is given by 

c = ^3 Jd 3 P 2 dn M (i-f^)(i + fi B) )fPfi B) 

-f (F) A B \l-fi F) )(l + fi B) )}, (4) 

where /« = /«(r,p,i) and fi j) = f^(r, Pi ,t), dQ f is 
the element of solid angle for the outgoing relative mo- 
mentum P3 — P4, v = |v — V2I is the relative velocity of 
the incoming particles, and a — Aira 2 is the scattering 
cross-section. The collision satisfies conservation of mo- 
mentum (p + P2 = P3 + P4) and energy (e + e 2 = £3 + ^4), 
with ej =p 2 j /2m j + TJV). 

The solution of Eqs. is carried out by using a nu- 
merical approach based on particle-in-cell plus Monte 
Carlo sampling techniques, which allows us to evaluate 
the kinetics of such systems down to T ~ 0.1 Tp. The 
technical details of the method have been given elsewhere 

In the following we shall focus on the dependence of 
the dynamics of a mixture of Np = 10 4 40 K atoms 
and Nb — 10 2 87 Rb atoms on the anisotropy A of the 
traps at fixed values of the average trap frequencies 
Cup = (up ± Ljp tZ ) 1 / 3 and ujb — i w b : z) 1 '' 3 - These 
are taken as the geometric averages of the trap fre- 
quencies in the experiments carried out at LENS on 
40 K- 87 Rb mixtures J2.2f, i.e. we set u> F = 2tt x 134.1 



s 1 and ujb = 2-7T x 91.2 s 1 . We fix the tempera- 
ture at T = 0.2 Tp, with Tp = Uuip {QNp) 1 ^ /k B be- 
ing the Fermi degeneracy temperature for noninteracting 
fermions. This also corresponds to T ~ 2.7Tbec where 
Tbec = 0-94:Tiujb {Nb) 1 / 3 /k B is the condensation tem- 
perature for the noninteracting Bose component. Finally, 
we assume a repulsive s-wave scattering length a = 2000 
Bohr radii. This choice of the fermion-boson scatter- 
ing length is dictated by computational convenience, but 
still leaves the mixture at A = 1 in a collisionless-to- 
intermediate scattering regime as for the 40 K- 87 Rb mix- 
tures studied at LENS. 



III. OSCILLATION MODES 

A. Dipolar oscillations 

The collisional state of the gas is revealed by the be- 
havior of the frequencies and damping rates of collec- 
tive excitations as functions of the collision frequency. 
Briefly, for dipolar modes the hydrodynamic behavior is 
signalled by a common frequency of oscillation of the two 
species and by a decrease of the damping rate with in- 
creasing collision frequency. Conversely the collisionless 
regime, which is attained at low collisionality, is charac- 
terized by different oscillation frequencies and by grow- 
ing damping rates. The collision rate can be evaluated 
either from a numerical simulation which actually counts 
the number of collisions at each time step, or by direct 
integration of the collision integral over momenta. Pauli 
blocking strongly quenches collisions at the temperatures 
of present interest and its handling in numerical studies 
requires suitably adapted methods of Monte Carlo sam- 
pling BIHEl- 

In the numerical simulation we excite dipolar modes by 
initially shifting the fermionic density profile by = 
yh] {mpOJF.i_) in either the axial or radial direction 
(z and x, say) and then fit the time evolution of the 
fermionic center-of-mass coordinates with the functions 
cos(fti t + <fi) exp(— ji t). From these fits we extract the 
oscillation frequencies Qi and the damping rates 7$ along 
the two directions. The results obtained from the numer- 
ical solution of the Vlasov-Landau equations are shown 
in Fig. ^ The error bars in the plot have been estimated 
from the standard deviations of fii and 7i in the fitting 
process. 

Starting from A = 1 and increasing the anisotropy to- 
wards an elongated cigar-shaped trap, the frequencies of 
the dipolar oscillations show different behaviors along the 
axial and the radial direction. The value of il±/u>p t ±_ in- 
creases slightly towards unity as A is decreased, whereas 
£l z /u)p jZ decreases with A towards an appreciably lower 
value. At the same time both damping rates tend to 
vanish, although that for axial oscillations appears to go 
through a broad maximum before doing so. These be- 
haviors indicate that in strongly elongated traps the dy- 
namics of the Fermi gas is collisionless in radial dipolar 
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oscillations, but collisional in axial ones. 

The coexistence of collisionless and hydrodynamic be- 
haviors in the dipolar oscillations as presented above is 
supported by the solution of scaling equations in the clas- 
sical limit. Within a classical model one can write a set of 
equations for the coupled motions of the centers of mass 
of fermions and bosons in the small-oscillation regime 
[l8L |22||. If we include mean- field effects, the equations 
of motion in the i-direction for the center-of-mass coor- 
dinates and for the relative coordinates are 

J icM,i = -^hd, 4 x CM,i ~ M r Auj 2 X r ,i/M 

\ X r<i = - AuJ 2 X C M,t - W, 2 j X r ,i + W 2 lf i X r ,i ~ X r>i 

(5) 

where icm.i = {vtifNfXf^ + mBNBXB,i)/M and x r ^ = 
XF,i — XB,i- We have defined Auj 2 — u> Fi — w Bi and 
u) 2 . t = {■mBN B UF,i + mFNFU 2 B i )/M with M = m F N F + 
tueNb and M r = msmp-A^sA^^/M. The hydrodynamic 
frequencies are given by 



motions we find 



m F N F uj f ,i + m B N B u 2 B i \ 
M 



(6) 



and the collisional frequency is luq = 4Q(itl b /Nf + 
vriF /Nb) I '(SrriF + 3ms), Q being the total number of 
collisions per unit time defined as 



W ^ k B T 



k x ky k z 



1/2 



(7) 



i) are the effec- 
rrij u) 2 j. In the 



In Eq. @ ki = k F ,ik B ,i/ '{k F ,i + k 
tive oscillator constants, with kj t j = 
large-Q limit Eqs. JSJ predict that the relative motion 
of the two clouds is overdamped, so that fermions and 
bosons oscillate together at the hydrodynamic frequen- 
cies. For the system parameters that we are using these 
are Whd,i = 0.9942 u> While due to the low number 
of impurities the value of cohd,i is very close to the bare 
trap frequency u>F,i> their difference can be amplified by, 
e.g., increasing the mass of the impurities or varying the 
trapping frequency. 

The mean-field correction to the bare frequency of the 
relative motion in Eq. (J5J is 



M r 



8n (F) 8n (B) 
dxi dxi 



(8) 



where 



AF,B) 



(r) are the equilibrium density profiles. 



Equation © extends the result of Ref. |23j to a general 
gaseous mixture. According to Eqs. JSJ the mean-field 
correction is more important in the collisionless limit and 
does not affect the value of the hydrodynamic frequency. 
However, it shifts the value of Q lock at which the two 
clouds become glued together. The locking point can be 
estimated by looking for an overdamped oscillation in the 
solution for the relative motion in Eqs. JSJ. If we neglect 
the coupling between the center-of-mass and the relative 
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(9) 



We have estimated a value of u)^ afi /uj Fi ~ 8 x 10~ 3 for 
the isotropic trap and verified that this ratio is slightly 
decreasing with A. Therefore the mean-field correction 
is negligible in the present case and the locking point is 
fixed by the bare frequency of the relative coordinate, 
namely WQ| Qlock ~ 2u r ^ ~ lAw F>i . 

The classical model does not contain the effects of Pauli 
blocking at low temperature, so that we have taken Q as 
a single fitting parameter at all values of A. The results 
obtained from Eqs. (JSJ with Q = 55ujf are reported in 
Fig. n an d compared with the numerical results for the 
dependence of the mode frequencies and damping rates 
on the anisotropy parameter. At A = 0.05 the gas is 
already very close to the hydrodynamic regime in its mo- 
tions along z, since the value of the oscillation frequency 
obtained in the simulation is close to u>hd,i in spite of the 
low concentration of impurities. The overall agreement 
between the simulation and the model is fairly good, and 
the deviations may be attributed to low-temperature ef- 
fects and to shape deformations of the distributions not 
entering Eqs. ©. 

Finally, it is worthwhile stressing that even though the 
effects on the dipolar oscillations due to the scattering 
with impurities are small, they can be greatly increased 
by tuning the experimental parameters. For instance, if 
we consider the parameters of the recent experiment by 
Takasu et al. |24j and take 100 174 Yb atoms as impu- 
rities, the difference between the oscillation frequencies 
along the axial and radial directions as shown in Fig. ^ 
will be of about 75% in the limit of small A. Therefore, 
we may expect that the effects here described will be- 
come easily observable in the near future when mixtures 
of 40 K and 174 Yb will be experimentally realized. 



B. Monopolar and quadrupolar oscillations 

Dipolar modes are a simple example of collective modes 
that can be experimentally analyzed and whose features 
expose the collisional state of the gas. As illustrated by 
the classical model, these modes can be viewed mainly 
as coupled oscillations of the centers of mass of the 
two clouds in the absence of appreciable deformations 
of their shapes. Other collective modes can be explored 
by deforming the clouds in different ways. The lowest- 
frequency modes of this type are the surface monopolar 
(I = 0) and quadrupolar (I = 2) modes, and these can 
be excited by small deformations of suitable symmetry. 
While in an isotropic trap the monopolar mode and the 
^2 = quadrupolar mode can be independently excited, 
in an axially symmetric trap the angular momentum is 
not a good quantum number and the two oscillations with 
l z = are coupled to each other. 
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FIG. 1: Fermion dipolar frequencies Qi (top panel) and damp- 
ing rates 7, (bottom panel) in units of WF,i, as functions of 
the trap anisotropy A. The results of the simulation (sym- 
bols with error bars) are compared with those obtained from 
Eqs. 0: open circles and solid lines denote radial oscillations, 
while full circles and dashed lines refer to motions along the 
z axis. The arrow in the top panel marks the locking point 
given in Eq. JSJ. 



As the collisionality of the fermion gas is increased, the 
transition to the hydrodynamic regime in strongly elon- 
gated traps manifests itself through changes in the £ z — 
mode frequencies from 2lqf,z to fihd = y^2/5u) FtZ an d 
from 2lu f , ± to y/lO/3u F , ± jUSEE]]]. A similar tran- 
sition can be expected in the present system, where the 
impurities act to transfer momentum between noninter- 
acting fermions. We have monitored several dynamical 
averages of the fermion cloud in the course of the simu- 
lation and in particular we have analyzed the time evo- 
lution of x(t) — (2v Fz — v F± ) after compression of the 
cloud density by about 5%. This average measures the 
anisotropy of the velocity distribution and is the respon- 
sible for the coupling between monopole and quadrupole 
excitations |27| . At variance from the dipole modes, in 
the monopole and quadrupole modes the radial and axial 
motions are strongly coupled and this requires that we si- 
multaneously follow the dynamics of the gas on two quite 
different time scales. In Fig. [21 we plot the value of the 
lower £ z = mode frequency obtained from the Fourier 
transform of x(t)- At very large anisotropy this peak fre- 
quency suddenly drops from the collisionless value 2u> F z 
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FIG. 2: Frequency Slo (in units of ujf, z) of the lower oscillation 
mode with t z — 0, as a function of A (symbols with error 
bars giving standard deviations). The solid line shows the 
lower i z — frequency for a interacting fermion system with 
&FT mq — 15 (see text). 



towards the hydrodynamic value a/12/5 wf ;2 , whereas 
the higher peak frequency (not shown in Fig. [2J stops 
at its collisionless value 2o>f,_l- Again collisionless and 
hydrodynamic behaviors coexist in a strongly elongated 
trap. 

The behavior of monopolar and quadrupolar modes 
with finite relaxation time T mq in a gas of interacting 
fermions inside an anisotropic trap has been studied by 
Vichi 25], who derived an implicit polynomial equation 
predicting a very steep downturn of the lower mode fre- 
quency in systems with ui F r mq S> 1 . The solid line in 
Fig. |21 shows the result of fitting Vichi's model to our 
data with the choice tu F T mq = 15. The sharp transition 
to the collisional regime is due to the fact that the impu- 
rities mediate the fermion-fermion scattering and lead to 
an effective relaxation time which is larger than that in- 
volved in the impurity-fermion scattering. Even though 
the effect is limited to large anisotropics, it could be- 
come observable at moderate anisotropics by increasing 
the numbers of particles or the strength of the boson- 
fermion repulsion. 



IV. EXPANSION DYNAMICS 

Many experiments extract information on the prop- 
erties of an ultracold trapped gas after it been allowed 
to undergo free expansion. This improves the spatial 
resolution of in- situ measurements and also gives access 
to the momentum distribution of the atoms. A fermion 
gas in the collisionless regime becomes spatially isotropic 
as it expands, regardless of the initial anisotropy of its 
density profile, whereas in the hydrodynamic regime the 
density profile inverts its aspect ratio during expansion. 
Furthermore, since the number of collisions diminishes as 
the gas expands, the hydrodynamic picture will become 
invalid when there are not enough collisions to sustain 
local equilibrium. 
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FIG. 3: Asymptotic aspect ratio R±/R z of the expanded 
fermion cloud as a function of the trap anisotropy A. The 
error bars are estimated from standard deviations in the fit- 
ting process. The solid line is the solution of Eq. I1H with 
to = 6/ 'Of. The collisionless and hydrodynamic limits are 
shown by the dashed and short-dashed lines, respectively. 

In the numerical simulation of the expansion, start- 
ing from the equilibrium density profiles of the trapped 
gaseous mixture we switch off the confinement and allow 
evolution according to the VLE. During the evolution we 
adaptively change the size of the computational domain 
in order to ensure that all particles are inside it. This per- 
mits us to study the free expansion of the cloud for long 
periods of time up to expansion ratios b{ (t) = Ri (t) /Ri (0) 
of a few hundreds, with Ri(t) — {x 2 F ^t)) 1 / 2 being the 
width of the cloud in the i-direction. In Fig. [3] we 
show the asymptotic value R±/R z of the aspect ratio 
of the fermion cloud as a function of the anisotropy 
A. The aspect ratio in the case of isotropic or strongly 
anisotropic confinement is equal to unity as in the col- 
lisionless regime, while for intermediate values of A the 
scattering against the impurities makes R±/R z deviate 
from unity. This behavior can be seen as a consequence 
of the competition between the time scales associated to 
axial and radial motions, as explained below. 

We have evaluated the total numbers Q± and Q z of 
collisions occurring during a lapse of time from t = to 
t± = 1/u>f,±. and to t z = l/u)p >z , respectively. Assuming 
that the collision rate Q(t) scales in time with the volume 
occupied by the gas and that the expansion dynamics is 
close to that of a collisionless system, Q± and Q z can be 
written as 

rtx,U 

Q±, z = / dtQ{t) (10) 
Jo 

= Q(0) / " ' dtil + ^^Y^l+ulj 2 )- 1 ' 2 . 
Jo 

In Eq. UTJ) we have taken Q(t) = Q(0)/(b 2 ± b z ) and set 
h(t) = (1 + ujj,^ 2 ) 1 / 2 28]. The integral in Eq. (TTJ 
is straightforward and yields the results shown in Fig. 
0] The following points should be noted: (i) the char- 
acteristic time t± for the radial expansion and conse- 
quently Qj_ decrease with increasing the anisotropy; and 
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FIG. 4: Numbers Q±/N F (solid line) and Q Z /N F (dashed 
line) of collisions per fermion, as functions of the anisotropy 
parameter A. 

(ii) the integrand of Eq. for Q z becomes negligible 
for tj_ <c t < t z as the fermion density drops due to rapid 
expansion in the radial direction, so that Q z increases at 
first with decreasing A, reaches its maximum at A ~ 0.4, 
and then rapidly drops as the density of the expanding 
cloud rapidly vanishes. 

The non-monotonic behavior of R±/R z as a function 
of A in Fig. [3] can be understood from the features of Q± 
and Q z in Fig. 21 The number Q± of collisions dimin- 
ishes with decreasing A while Q z increases, so that col- 
lisions play different roles in the axial and radial expan- 
sion. Both Q± and Q z are dropping at large anisotropics, 
so that collisionless behavior is emerging in the ballistic 
expansion from a strongly elongated trap. 

The behavior of the expanding g elS cLS el function of A 
can also be analyzed by means of scaling equations. In 
this approach |28| the expansion of a fermion gas with 
finite collisionality can be described by the equations 

( 2 0i 

where 6i (t) are effective temperatures along the two spa- 
tial directions, T{9 i ,b l ) = r (fX,- J2k #fc)~ 1/2 and T o 
is the collision time at t — 0. The collisionlcss and hydro- 
dynamic limits correspond to taking r = oo and r = in 
Eqs. (flip, respectively. 

We have taken the collision time r in Eqs. i|llfl as a fit- 
ting parameter to represent the indirect fermion-fermion 
scattering induced by the impurities. The solution of 
Eqs. with the choice To = 6/lof reproduces the non- 
monotonic behavior of the aspect ratio as shown in Fig. 
13 Finally, we have checked that the qualitative behav- 
ior shown in Fig. ^for Q± and Q z does not depend on 
the assumption of a collisionless evolution of the gas and 
persists on describing its expansion with any value of t . 
The magnitude of the deviation from unity of the aspect 
ratio of the expanding fermions depends on the value of 
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To and can be increased by changing the fermion-impurity 
interaction via Feshbach resonances | 29j or by perform- 
ing^ an out-of-equilibrium experiment as described in Ref. 

V. SUMMARY AND CONCLUDING REMARKS 

We have studied the low-lying surface modes and the 
ballistic expansion of a spin-polarized Fermi gas interact- 
ing with thermal impurities as functions of the anisotropy 
of its confinement in a cigar-shaped harmonic trap. We 
have solved for this purpose the Vlasov-Landau equa- 
tions for the dynamics of the mixture and compared the 
results with simple scaling equations containing collision- 
time fitting parameters. The results show that for large 
anisotropies a collisionless behavior in the radial dipo- 
lar oscillations and a hydrodynamic behavior in the axial 
ones are simultaneously established. For monopolar and 
quadrupolar excitations we have observed a collisionless 
spectrum irrespectively of the strength of the anisotropy 
parameter, except for extremely large anisotropies where 
the frequency of the lower l z = mode decreases towards 
the hydrodynamic value. 

On the other hand, during ballistic expansion the two 



different time scales for radial and axial motions enter 
into competition. The result is that the expansion of a 
strongly anisotropic cloud is essentially collisionless due 
to the rapid drop of the particle density, whereas at in- 
termediate values of the anisotropy the aspect ratio is 
sensitive to the collisions. 

The analysis presented here has only concerned slightly 
doped fermions in a collisionless-to-intermediate scatter- 
ing regime as, e.g., that attained in the experiments at 
LENS. The effects of the trap anisotropy could be fur- 
ther enhanced exploiting the rich variety of experimental 
set-ups and trapped isotopes available. It would be also 
interesting to extend our work to strong-interaction situ- 
ations among fermions and impurities or among fermions 
in two-component Fermi gases approaching the unitary 
limit. 
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